rm(list = ls(all = TRUE)) #clear workspace

library(foreign)
library(survey)
options(survey.lonely.psu = "adjust")

setwd("C:/Users/martin.romero/Dropbox/ENVIPE/2009")

################################################
#Pregunta al entrevistado
################################################

data <- read.dbf("delitos_2.dbf")

#Ponderador y diseno
data.design <- read.dbf("delitos_3.dbf")
data.design <- data.design[,c("FAC_ELE_B","EST","UPM")]
data <- cbind(data, data.design)

data$FAC_ELE <- as.numeric(as.character(data$FAC_ELE_B))
design.per <- svydesign(id=~UPM,strata=~EST, weights=~FAC_ELE, data=data)

#4.1 Mencion� que en 2008 usted fue v�ctima de alg�n delito. 
#�Podr�a decirme de cu�l o cu�les delitos se trat�?
#Secuestro o secuestro expr�s
table(data$DP4_1_4, useNA="ifany") #34 casos
data$victima <- factor(0, levels=c(0,1))
data$victima[data$DP4_1_4=="X"] <- 1
table(data$victima) #check
design.per <- svydesign(id=~UPM,strata=~EST, weights=~FAC_ELE, data=data) # ...again
foo <- svytotal(~victima, design=design.per, na.rm=TRUE) #51,049
print(foo)
confint(foo) #24,918-77,179
cv(foo)

#4.2 �Cu�ntas veces le sucedi�?   					
table(as.numeric(as.character(data$DP4_2_4)))
sum(as.numeric(as.character(data$DP4_2_4)), na.rm=TRUE) #34
foo <-svytotal(~as.numeric(as.character(DP4_2_4)), design=design.per, na.rm=TRUE) #51,049
print(foo)
confint(foo) #24,918-77,179
cv(foo)
